A Random Matrix Models of Relaxation 



m 
o 
o 



J. L. Lebowitz 
Department of Mathematics, Rutgers University, USA * 

L. Pastur 

Department of Mathematics, University Paris 7, France ^ 



Abstract 



We consider a two level system, S2, coupled to a general n level system, Sn, via a 
CN I random matrix. We derive an integral representation for the mean reduced density 

matrix p{t) of ^2 in the limit n — > oo, and we identify a model of 5„ which possesses 
^ ■ some of the properties expected for macroscopic thermal reservoirs. In particular, 

! it yields the Gibbs form for p(oo). We consider also an analog of the van Hove limit 

^ I and obtain a master equation (Markov dynamics) for the evolution of p{t) on an 

■ appropriate time scale. 
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1 Introduction 

^ ' The time evolution of a "small" quantum system interacting with a "big" thermal reservoir 

is a much studied problem in various contexts, see review works [3 13 01 IHl El El Cj ■ Two 
natural questions for such models are the large time limit of the state of the small system 
^ . and the nature of the approach to this limit. It is generally expected and rigorously proven 

for certain types of reservoirs and couplings that, starting from an arbitrary initial state 
of the system and an equilibrium Gibbs state at temperature T of the reservoir, the final, 
t ^ CO, state of the system is given by a projection of the joint system plus reservoir 
Gibbs measure at the temperature T[Tl|Uj. 

In this paper we discuss the above questions in the frameworks of random matrix the- 
ory. Namely, we consider a small system S2, having only two levels (like in the archetype 
spin-boson model [T71 HJ IHj), coupled to a system with n levels (denoted by 5„) via an 
n X n random matrix. Similar models have been discussed in papers [21E1CSIE1- Here 
we will investigate in a more explicit way what features of the n-level system will lead, in 
the limit n 00, to the reduced density matrix p{t) of the small system having a Gibbs 
form for large times and whether p{t) is described by a Markov evolution. Our analysis 
will be based entirely on the fact that the interaction is given by a typical realization of 
a random matrix. This models in some way a strongly interacting reservoir. 

The paper is organized as follows. In Section 2 we describe the model and define 
the reduced density matrix for the small, two level, system. Since it is random, its full 
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statistical description is of interest. In this paper we confine ourselves to the analysis of its 
mean which we argue is also its typical behavior in certain cases. In Section 3 we consider 
the equilibrium properties of the combined system. We find the reduced density matrices 
of 1S2 corresponding to the microcanonical and canonical description of the composite 
system and give conditions for their equivalence. 

Section 4 deals with time dependent (non-equilibrium) properties of the model. We 
find that the t —>■ 00 limit of p{t) depends in general on the initial density matrix. However, 
for the reservoir introduced in Section 3 designed to guarantee the equivalence of the 
microcanonical and canonical description in equilibrium, the t ^ 00 limit of the diagonal 
entries of the mean reduced density matrix is independent of the initial density matrix 
and has the Gibbs form. In Section 5 we treat the van Hove asymptotic regime (small 
coupling, long times). Here both the equilibrium and the time dependent form of the mean 
reduced density matrix can be expressed directly via the density of states of the reservoir. 
We find that for the special form of the density of states of iS„ of Section 3 the t — ^ 00 
limit of p{t) is the Gibbs distribution, and its time evolution is Markovian. However, in 
this case the obtained time decay rates prove to be vanishing. In Section 6 we consider 
a more general model, where one can obtain quite reasonable decay rates. The model is 
similar to that, studied in [T2]. In Section 7 we continue studies of the time evolution p{t) 
discussing its dependence on the initial density matrix and on the parameters. Several 
formulas, used in the main body of the paper, are derived in Appendices. 



2 Description of Model 

Let hn be a real symmetric nxn matrix with eigenvalues Ej, j = 1, n. We characterize 
the spectrum of /i„ by its density of states 

n 

ut\E)=n-^Y.^{E-E,). (2.1) 

i=i 

We assume that i/q"'' converges as n ^ cxo to a limiting density vq, i.e. that for any 
continuous and rapidly decaying function ip we have: 

/oo 1*00 1*00 

ip{E)v^^\E)dE = (p{E)iyo{E)dE, / iyo{E)dE = 1. (2.2) 
-00 J —00 J —00 

Let Wn be a real symmetric nxn random matrix, whose probability density is 



g„^exp |-^Tr w^^ 



(2.3) 



where Qn is the normalization constant. In other words, the entries Wjk, ^ < j < k < n 
of the matrix w„ are independent Gaussian random variables, with 

=0, {w%) = {l + 5,,), (2.4) 

where the symbol (...) denotes averaging with respect to the distribution ()2.3p . This 
probability distribution is known as the Gaussian Orthogonal Ensemble (GOE) [T3|. 
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We define the Hamiltonian of our composite system S2^n as a random 2n x 2n matrix 
of the form 

Hn = sa' ® 1„ + I2 ® /in + va"" ® Wn/n^^'^, (2.5) 
where 1; {I = 2, n) is the I x / unit matrix, a^, and are the Pauli matrices 

The symbol denotes the tensor product, and s and v are positive constants. 

The first term in ()2.5p is the Hamiltonian of ^2, the second term is the Hamiltonian 
of Sn, and the third term is an interaction between them. Thus s determines the energy 
scale of the isolated small system, and v plays the role of the coupling constant between 
the two level system S2 and the n level system iS„. Note also the n~^/^ scaling of the 
interaction. 

We emphasize that letting n — 00 does not imply that iS„ will behave like a macro- 
scopic system. iS„ could perfectly well be just a single particle moving on a periodic chain 
of n sites with spacing 6, for which, in the limit n 00, 

jjo{E) = (n^E{A6-^ -E)^ \ 0<E < 45-\ (2.6) 

In fact the normalization assumption ()2.2jl implies that we are not dealing here with what 
one usually thinks of as a thermal reservoir. 

We are going to study the time evolution of the mean reduced density matrix of ^2, 
assuming that at t = the density matrix of the composite system iS2,n is of the form 

;,W(o) = p(o)®P,, (2.7) 

where p(0) is a 2 x 2 positive definite matrix of unit trace and Pk is the projection on the 
state of energy Ek of the reservoir. Let 

/^^"^W = {/^SUW, a,/3 = ±, j,fc = l,...,n} (2.8) 

be the density matrix of the composite system S2,n at time t, corresponding to the initial 
density matrix /im^(O) of ()2.7|) . Then the mean reduced density matrix of the small system 
is defined as 

n n 

&E„t) = J2 {f^'aUt)) = E ((^'"''>^^"H0)e^*^")„^.,J . (2.9) 

We have evidently for all k = 1, ...,n and t > : 

E/^S(^'^'^) = 1- (2-10) 

a=± 

The mean reduced density matrix is determined by the initial density matrix p(0) of S2 
via the linear relation 

^Ek^t) = J2 tts(Ek,t)p,,,{0), (2.11) 
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where 

n 

i=i 

is the "transfer" matrix. It can be viewed as an analogue of the influence functional, 
introduced by Feynman and Vernon in a similar context [H]. 

A particular case of the mean reduced density matrix is the mean transition probability 
p^^l{E,t) for ^2 to be found in the state |a) at time t provided that at t = it was in 
the state \P) and the initial state of Sn was the pure state of energy for some k. It is 
given by 

P^:iiE.,t) = Til(E,,t) = X: (|(e-*^")^^.,,f ) , (2.13) 

j=i \ I 

with 

E^'S(^^'^) = 1' '^ = ±- (2-14) 

a=± 

The diagonal entries p'^)y{Ek,t), a = ± are known as survival (or return) probabilities. 

We can evidently replace Pk by any superposition Ylik'^kPk of pure states of iS„, e.g. 
by that, corresponding to the canonical density matrix e"'^'^" /-^n, where Zn is the par- 
tition function of iS„ and (3 is its inverse temperature. The initial density matrix for the 
composite system will now be (cf (j2.7j) ) 

/^i"H0) = MO) ® e-^'^" /Z„ , (2.15) 

corresponding to the choice 

Ck = e-^^^ In r e-^^z/J")(E)dE , (2.16) 

/ J — oo 

where is defined in fl2.1|) . 

It should be mentioned, that our results remain valid also in the case, where the 
matrix /i„ is non-diagonal, random, independent of Wn, and possesses property ()2.2p with 
probability 1. However, the randomness of plays practically no role. Besides, our 
results are not limited to the special Gaussian distribution (j2.3|) - (j2.4|) of the interaction, 
but hold also for any real symmetric random matrix Wn in ()2.5|) . whose entries Wjk, 1 < 
j < k < n are independent, verify ()2.4|) and have fourth moments bounded uniformly in 
n. The techniques in this case are, however, more involved. 



3 Equilibrium States 

3.1. General formulas. 

We begin by considering the equilibrium (time independent) microcanonical density ma- 
trix of the composite system S2,n- 

n{X) = 6{X - Hn) /Tr 6{X - H^) , (3.1) 

Following a standard prescription of statistical mechanics, we will replace the Dirac delta- 
function in ()3.1|) by the function {2e)~'^Xe, where Xe is the indicator of the interval {—e, e), 
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and e <^ X. Then the reduced micro canonical density matrix of S2, i.e., the microcanonical 
density matrix of S2,n, traced with respect to the states of Sn, is the 2x2 matrix of the 
form 

^-(A)^^^, (3.) 

where 

n 

z7W(A) = (2en)-i ^Xe(A - H^U,,. (3.3) 
i=i 

The corresponding canonical distribution of the composite system is 

and the reduced distribution of the small system is 

, (3.5) 

E..±ne-/^V(})(A)dA 

where (cf dSH) 

v^::i{\)=n-^Y.^{\-H^U,^. (3.6) 
i=i 

By using standard techniques of random matrix theory (see e.g. [HI El), it can be shown 
that with probability 1 with respect to the distribution ()2.3|1 of the interaction matrix the 
n — > 00 limit of ()3.3j) is 

= Sa,^Va, z7q(A) = (28)'^ / l'a{.lj)dll (3.7) 

with 

i/^(A) = 7r"^$5r„(A + iO). (3.8) 
The pair ra{z) solves the system of two coupled functional equations 

ra{z) = / 2 rr, a = ±, 3.9) 

with uq, defined by ()2.H) - ()2.2|) . The solution is unique in the class of functions analytic 
for Qz 7^ z, and satisfying the condition '^ra{z) ■ Qz > 0, z 0. 

The n ^ oo limit of the reduced microcanonical distribution (|3.2p is then 

Cl;„^(A) = Sa,^uJaW, UJa{X) = _ ^ ^■."'['^l (3.10) 

Correspondingly, the n — > oo limit of the reduced canonical distribution ()3.5|) is 

r e-^V„(A)rfA 

The absence of off-diagonal elements in p{j3) is due to the special structure of the inter- 
action (proportional to the matrix). This is a rather common phenomenon in random 
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matrix theory and many ensembles, studied so far, have this property. If, on the other 
hand, one puts an arbitrary 2 by 2 matrix instead of cr^ then p{f3) will be non-diagonal. 

In the next section we will consider models of iS„ in which ()3.10p and ()3.1H) are the 
same for the proper choice of (3. It will turn out that we will have not only the "equivalence 
of ensembles" but also that the Pa{f3) will be independent of the interactions and therefore 
be of the Gibbs form, 



2 cosh Ps 



(3.12) 



3.2. A model for the reservoir. 

We will consider now a model in which Sn consists of J particles or other subsystems, 
each having m levels (orbitals). We then have n = m"^. We assume further that the 
subsystems do not interact. The (normalized) density of states z/q"^ (12.11) is then the J- 
fold convolution of the (normalized) density of states of a subsystem, and we assume 
that in the limit m — ^ oo will approach some normalized density q. 

Following the limit m —>■ oo we will take the limit J —>■ oo. Note that this is very 
different from the usual thermodynamic limit, e.g. for J particles on a lattice of m sites, 
where J would grow with m in such a way that we would have for large J that the number 
of levels between E = Je and E = Je + dE grows as exp{JSj{e)}dE. The entropy Sj{e) 
would then be a monotone increasing function of e which would approach, as J oo, 
the entropy/per particle, ^(e) with S'{e) = /3, the reciprocal temperature of the reservoir. 
The same happens in spin systems where m remains finite as J — *■ oo, in which case S{e) 
would have a maximum at some finite e, after which it would decrease. 

The reason we cannot follow this procedure here is the dissonance between random ma- 
trix theory and thermodynamic behavior already mentioned above. To get a (normalized) 
density of states for the thermodynamic system would require that exp{JSj{E/ J)}/m'^ 
have a limit as J ^ oo which is not the case. 

We are therefore forced to keep J fixed while taking the limit m ^ oo, yielding a 
continuous q and z/q as its J-fold convolution. This form of z/q is standard in probabilities 
rather than in statistical mechanics. However, they both suggest that we consider cases 
where 

lim sj{e) = s{e), sj{e) = J^Mog z/o( Je). (3.13) 

J^oo 

The function s, appearing above, is in fact the rate function of the large deviation theory 
(see e.g. jl]) and can also be seen as an analog of the entropy of statistical mechanics, 
although, since is normalized to unity, the analogy can not be too close. In particular, 
s will be negative everywhere, since ^^{Je) ^ as J ^ oo, unlike the genuine entropy in 
quantum statistical mechanics. 

Formula (|3.13p suggests the following one: 

z/o(Je + e)/vQ{Je) = exp{ J[sj(e + e/J) - sj{e)]} e''(")^ J ^ oo. (3.14) 

Hence, the equahty 

s'{e)=P, (3.15) 

analogous to the well known statistical mechanics relation between the entropy and the 
temperature, implies the formula 

i/o(Je + e)/i/o(Je) ^ (3.16) 
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To obtain a positive (3 we have to choose an interval on which s(e) of (|3.13|) increases. 

However, for the analysis of (|3.1Up . and (|3.11|) . as well as in the study of the time 
dependent properties of the model, we will need the asymptotic formula 

z/„(Je)/i/o(Je) ^ e"^"", J ^ oo, (3.17) 

where e and 13 are related by ()3.15|) . This requires certain conditions on q (in particular, 
the existence of certain exponential moment) and a certain amount of technicalities. To 
keep our presentation simple we consider a model case, where q is the Gaussian density 

^^^^ = (27raV/^ ^^P {-(e - eo)V2a'} • (3-18) 
In this case we have immediately 

""'^^^ = (27rjl^)V2 - ^eo)V2Ja2} , (3.19) 

and, hence, in view of (j3.13j) . 

s(e) = -(e-eo)V2a^- (3-20) 

Furthermore, we prove in Appendix 1 that if is given by ()3.19|) then the density z/^, 
defined by ()3.8p - ()3.9p . verifies the asymptotic relation ()3.17p . 

It is easy now to show that this model yields a reduced density matrix of S2 that is 
diagonal in a and has the same Gibbs form (|3.12j) for both the microcanonical and the 
canonical cases: 



lim lim 



^^j^ J-00 m-00 ^^^^ e-^^V^"^ {Je)de 2 cosh (3s ' 

(3.21) 

provided that the energy e "per particle" and the temperature (3 of the reservoir are 
related via the usual equation ()3.15p . The first formula follows immediately from ()3.17p . 
To obtain the second formula we have to take into account that fj3.15|) is the saddle point 
condition for the respective integrals in the n = 00 version ()3.11|1 of the canonical reduced 
density matrix ()3.5j) . 

We see that because of the rather "non-macroscopic" (in the thermodynamic sense) 
nature of the density of states the random matrix model of the reduced equilibrium 
description does not exhibit, in general, properties that are natural to expect from the 
statistical mechanics point of view. The same is valid for the time evolution, as will be 
seen in the next sections. However, the above ansatz for the density of states (be it given 
by (j3.19p or something different) leads to the equivalence of the reduced microcanonical 
and canonical description in the J — 00 limit of our composite system, if we assume, in 
addition, a natural from the statistical mechanics point of view relation ()3.15|) between 
the energy and the temperature. 

4 Time Evolution 

4-1. General formulas. 
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Assume that the energy E^. in (j2.12|) belongs to a subsequence {E^^} that converges to 
a given energy E of the reservoir such that vo{E) > 0. By using an extension of the 
techniques, presented in |0 [T^ , one obtains the following expression for the n —>■ oo limit 
of the transfer matrix ()2.12|] : 

T^p^s{E,t) = (4.1) 



(27rz) 




c Jc 



where the contour C encircles the real axis, 



ral3'y5{E,Zi,Z2) = lim ( R^k,5j{zi)Raj,l3k{z2) ) , (4.2) 



n—*oo 

\j=i 



and R{z) = {Hn — z)^^ is the resolvent of the Hamiltonian ()2.5|) . The "two-point" 
functions ra/3'ys{E, Zi, Z2) are analytic in zi and in Z2 outside the real axis and have the 
form 

raf3^siE,zi,Z2) = r^{E,zi)rp{E,Z2)D:^^p{zi,Z2) (4.3) 

X {6a,l3S^,5 + V^6-a,f36^^^sr-'y,-l3{zi , Z2)) 

Here the "one-point" functions ra{E,z), a = ± are 

^aiE, z) = —- i (4.4) 

E + sj — z — v^r_a\z) 

ra{z)^ a = ± solve system ()3.9|) . 

E>^,p{zi, Z2) = 1- v\^p{zi, Z2)r_^_fs{zi, Z2), (4.5) 

and 

r^A^iZ2) = J r^{E,z,)rp{E,Z2)ME)dE. (4.6) 

is another two-point function. 

These formulas, and 1)2.111) lead to the following expression for the n oo limit of 
the reduced density matrix 

PaAE,t) = j^dz2 j^dz.e''^^'-''^ (4.7) 

^ rgjE, Zi)rs{E, Z2)Pa,5(0) + V^r_a{E, Zi)r_s{E, Z2)ra,^(^l, Z2)p-a-5{0) 

1 - v'^ra,s{zi, Z2)r_a-5{zi, Z2) 
Hence the limiting transition probabilities, corresponding to Pa,a(0) = (5a/3, are 

Po.AE,t)= p^AE^t). (4.8) 

4-2. Reduced density matrix and the transition probabilities in the infinite time limit. 
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By using the analyticity of the integrand of (|4.7p in zi and in Z2, we can write the following 
integral representation for the diagonal entries of the reduced density matrix of the system 



oo 



p^AE,t) = -TTT-T^ / dX2 / rfAie**^^ (4.9) 

aP—a,—a (0) 

Here the integration path in Ai encircles A2 from below in the clockwise direction, 5 A = 
Ai — A2, 

5r^{E) = r^{E,\i-iO) -r^{E,\2 + iO), (4.10) 

and ra{z) and ra{E,z) are defined in ()3.9|1 and in ()4.4|) . 

These formulas imply that the diagonal entries of the reduced density matrix are, in 
the limit t 00, given by 

p„,„(E,oo)= f cu„(A)(^z/^(E,A)p^,^(0))rfA, (4.11) 

where 

^aiX) = ..'("f ^ a = ± (4.12) 
z/+(A) + z/_(A) 

is the e — form of microcanonical distribution of the small system, defined in fl3.1()|) 
with z7^(A) of (jSHI), replaced by z/^(A) of - (jHSl), and 

i/„(E,A) = -$5r„(E,A + iO) (4.13) 

TT 

t;V_„(A) 

— a = ±, 



(E + sa - A - t;2gfJr_„(A + iO))^ + 7rVz/2^(A) ^ 



where the last equality follows from ()4.4|) . The functions Va{E.,X) are non- negative and 
satisfy the relations 

/oo /»00 /"OO 

z/„(E,A)dA= / Uo.{E,\)dE =1, I u^iE, X)uo{E)dE = u^^X), a = ±. (4.14) 
-00 J —00 J —00 

Relations ()4.8|) and ()4.1ip imply that the transition probabilities are given by 

/oo 
uj^iX)u^{E,X)dX, (4.15) 
-00 

where Ua is the microcanonical equilibrium distribution ()4.12|1 . and i^aiE, A) is defined in 
()4.13|) - ()4.14|) . This expression satisfies the normalization condition ()2.14|) . 

As already remarked at the end of Section 3.1, there is no reason to expect that 
p{E, 00) will correspond to any equilibrium state for a general z/q. In fact we can show 
(see Section 7.1 ) that the diagonal entries pa^a{E, 00) of (j4.11|) will in general depend 
on p(0). The situation is different for the model of the reservoir, introduced in Section 
3.2. Indeed, setting E = Je{/3) in (jmH), where e{/3) is defined by dSlSI), and iTM . 
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we can change to the variable A = Je, and replace asymptotically the density Va{Je) by 
UQ^Je) exp{(e — eQ)sa/a?}, and the density UaiJ^iP), Je)J by 6{e + a? (3 — Cq) for J ^ oo, 
according to relations ()A.9|) . and ()A.10|) . As a result we obtain, in the infinite time limit, 
the canonical Gibbs distribution form ()3.2H) of the diagonal entries of the mean reduced 
density matrix of the small system with a temperature, determined by the energy by Sn- 
Note that the interaction does not appear in the limit although it ensures the loss of 
memory of p(0). 



5 Van- Hove Limit 

In this limit the coupling constant v of the system-reservoir interaction tends to zero, the 
time t tend to infinity while the transition rate, given by first order perturbation in the 
interaction, is kept fixed. In terms of ()4.H1 this corresponds to letting 

V ^ t oo, r = tv'^ fixed (5-1) 

after the limit n —>■ oo. 

By using general formulas (j4.7j) - (j4.3|) . it is easy to show that the off-diagonal entry 
p+_ of the reduced density matrix vanishes in the van Hove limit. 

Let us consider the diagonal entries pQ,^a{E,t) in this limit. Changing variables to 
A2 = A, 6\ = in the integral representation ()4.9p . and using ()4.10p . ()3.8p . and ()4.13p . 
we obtain 

'^^a(^)lAi=A2=A = -27rzz/„(E, A), Sr^l^^^^^^^ = -27iiu^{X), (5.2) 

where now 

z/^(A) = z/o(A - sa), iy^{E, A) = 5{E + sa - A). (5.3) 
The limit p"J^^{E^ r) of the diagonal entries, is then given by 



00 poo 



pll{E,r) = ^ / dX I di^^^ 



—00 J —00 



X 



faiE, A) - 27rzz/a(A)]pa,a(0) - 27rzz/_„(E, A)z/a(A)p_a _a(0) 



e(e-^r(A)) 

where 

r(A) = 27r(z/+(A) + z/_(A)), (5.4) 

and the integration with respect to ^ is taken by encircling the origin from below in the 
clockwise direction. Computing the integral with respect to ^ by residues, we get (cf 
Appendix 2) 

/oo 
U(A)^z/,(i?,A)p,,,(0) (5.5) 

^ z/»(^, A)z/_a(A)p»,a(0) - z/_a(E, A)^/a(A)p_^_^(0) ^_^^(A)^^^^ 

z/+(A) +z/_(A) ^ ) 

We see that the time independent part of the formula coincides with the r.h.s. of (j4.1H) . 
modulo the replacement ()5.3p . 
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Now, taking into account (|5.3p . we obtain finally that the reduced density matrix in 
the van Hove limit is 

pf,(E,r) = 5„,,pfjE,r), (5.6) 

where 

pS(i^, r) = ,^"'^11^, , P.AO) + ''r\^/ -.-.(0) (5.7) 

UolEj + UolE + Zsa) Uo[L - 2saj + UolE) 

^uo{E) + ME + 2sa)^"'"^''^^ ME - 2sa) + z/o(i5)^'"'""^ ^ 

and 

r„(E) = T{E + sa) = 271 {uo{E) + vo{E + 2sa)) , (5.8) 

with r(A), given by ()5.4p . In a general case of an arbitrary density of states uq of the 
reservoir the infinite time part of the reduced density matrix ()5.7p in the van Hove limit 
depends on the initial density matrix p(0) of the system, like it was for our general 
formula 1)4.111) (see, however. Section 7.2). Besides, the r.h.s. of formula ()5.7)1 contains 
three "modes": the time independent term, and the two terms with different exponentials 

in t. This implies that the pair (p++,p ) can not obey a system of two differential 

equations of the first order, an expected form of the master equation for p^^^{E, r) 

In fact, it is easy to see that if the first line of ()5.7|) is independent of p(0) for all values 
of E (recall that this is a free parameter of our model), then 

uo{E) = Ce^^, (5.9) 

where C and (5 are constants. This form of z/q leads evidently to the Gibbs form of the 
equilibrium distributions ()3.1(J|) . and to the right side of ()5.7j) having the form 



with 

Tc,{E) = A7iuo{E + sa) cosh Ps. (5.11) 

We note, however, that while the form ()5.9|) of uq for all E is incompatible with the 
normalization condition ()2.2|) . it can be valid "locally" in E as in formula ()3.16|) . We 
have seen in Section 3.2 that ()3.16|) implies the Gibbs form of the equilibrium densities 
matrices ()3.10|) . and p.lip in the zero coupling limit. 

Likewise, by using the asymptotic relation ()3.1fij) in the time-independent term of 
(15. 7j) . we obtain the first term of ()5.1())1 i.e. the Gibbs form of the time-independent term 
(15. 7p . after the replacement E by Je{P) and the subsequent limit J — > oo. 

The situation with the time dependent part of the reduced density matrix in the 
considered regime ()5.1|) - ()3.14|) is less satisfactory. Indeed, according to relation ()3.19p . 
the density of states z/o(Je(/?)) is exponentially small for large J (except the case, where 
e(/3) — Co is of the order J~^/^, that corresponds to the "infinite temperature" /3 = 0). 
Hence the decay rate (|5.11|) F will also be exponentially small in J. 

A simple way to make Ta{E) in ()5.8|) non- vanishing as J — > oo is to assume that the 
coupling constant v in ()2.5|) depends on J. This can be done if we write v as the product 
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of two factors Vi and f2, use Vi in the definition (|5.1|) of the van Hove hmit and keep V2 
in the decay rate. This yields 

V^{E) = AtivI{pq{E) + uo{E + sa)) (5.12) 

instead of ()5.8j) and allows us to rescale the decay rate by setting t>2 proportional to 
{vq{E) + vq{E + sa))~^. Since the parameter E appears in the initial formulas ()2.13|) and 
()2.12|) via the condition Ek = E, we see that we have to make the interaction term in 
()2.5p dependent on the energy levels of the reservoir. In the next section we consider a 
version of the model, possessing this property. 



6 A More General Interaction Model 



In this section we consider a model, given by the Hamiltonian ()2.5|1 . in which the system 
- reservoir interaction Wn is a real symmetric Gaussian random matrix with statistically 
dependent entries whose covariances are given by the relations (cf ()2.3j) ): 

E{wjfc} = 0, E{wjkWim} = f{Ej,Ek){Sji6km + SjmSki), j,k = l,...,n, (6.1) 

where f{E, E') is a non-negative symmetric function. Assuming that f{E, E') vanishes if 
\E — E'\ > b, for some b, we obtain the so-called band matrices, whose entries are non zero 
inside the "band" of the width b. The case f{E,E') = 1 corresponds to ()2.3p . Similar 
models were considered in [TH UHl 121 

A somewhat tedious calculation leads to a transfer matrix Tai3^s{E,t) that is given 
again by formula (|4.H) in which now (cf (j4.3j) ) 

rap^siE,Zi,Z2) = r^{E,Zi)rp{E,Z2) (6.2) 

X (^1 - S^a,-5{Z1, Z2)Sa,5{zi, 2:2)) 

X {5a,/35-y^S + 5-a,/35^^,5S-^-y-f3{zi, Z2)), 

where (cf (jOjl ) 

ra{E, z) = -— \ ,^ y (6.3) 

E + sa — z — A^a{E, z) 

The "self-energies" Aa{E, z), a = ± solve the system of two nonlinear integral equations 

(cf dniD) 

Sa,6{,zi, Z2) is an integral operator, acting on a function of by the formula 

{S^,s{zi, Z2)(t>) (E) = / v^fiE, E')r^{E', z^E', Z2)(l){E')v^{E')dE\ (6.5) 



and 

Sa,&{zi, Z2) = iSa,5iZi, Z2)1){E). (6.6) 

It is easy to see that if f{E, E') = 1, then the above formulas coincide with our previous 
formulas (lO) - in which A„(E, z) = v'^r^iz). 
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We show in Section 7.3 that particular asymptotic cases of the above formulas corre- 
spond to the results of papers jTH [12] . 

Consider now the van-Hove limit of this model as defined in ()5.1|) . One obtains then 
the expressions ()5.7|) - ()5.8|) in which z/q,(A) = z/o(A — sa) is replaced by 

z/o(A — sa)/(A + sa, A — sa). (6.7) 

Following ^21, we choose the form-factor / in ()6.1|) in the form 

/(A, /i) = = 

which is in agreement with our discussion at the end of the previous section. This leads 
to the following form of the diagonal entries of the reduced density matrix (cf ()5.1U|) ): 

n (F^\ e-/^- , e^^y^AO) - e-^-p-.,-.(0) ^_.r 
where now (cf fIS.lip ) 

r = 47r«;(2s)cosh/5s. (6.10) 

Unlike (|5.7|) . obtained in Section 5 for / = 1, the r.h.s. of (j6.9p contains only two 
"modes", the time independent term and the term with the exponential e~'^^. As a 
result, pa,a{E, r), a = ± is the solution of the system of two differential equations: 

^p,,„ = -2ttw{2s) (e^^>„,« - e~^^V_,,_,) . (6.11) 

Eq. (j6.1ip is a version of the master equation which is believed to describe, on an ap- 
propriate time scale [31^1 EI? the approach to equilibrium of a small system interacting 
weakly with a reservoir. An analog of 1)6.111) was obtained in |12| (see Section 6, in par- 
ticular formulas (52) - (53) there). From the probabilistic point of view ()6.11|) describes 
the two-state Markov process, known as the random telegraph signal or the dichotomy 
process (see e.g. jlD]); see discussion after (7.14). 



7 Some Properties of the Models 

In this section we discuss certain, mostly asymptotic, properties of the models, specified 
in Sections 2 and 6, that do not require the special structure of the reservoir, introduced 
in Section 3.2. 

7.1. The dependence of the reduced density matrix on the initial density matrix. 

We consider here the dependence of the t —>■ oo limit ()4.11|) of the diagonal entries of the 
reduced density matrix on the initial density matrix p(0) of the small system in the model 
of Section 2. 

It is easy to find from ()4.1ip that if the pa,a{E, oo) are independent of p(0), then they 
have the form 

1 r i^+(^,A) + z/_(E,A) 



Pa,a{E,oo) = - I — iya{X)dX, (7.1) 



oo 
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corresponding to p±(0) = 1/2 in (j4.1ip . and that if the pa,a{E, oo) are independent of 
p(0) for all E, then 

K(A)-..(A))!,,^„^ (7.2) 



. i.+(A) + ,._(A) 

Thus, the equality z^+(A) = z^-(A) has to be valid for all A which is impossible for the 
Hamiltonian ()2.5p with s ^ 0, i.e. if the small system has a nontrivial dynamics. Hence, 
there will in general exist energies of the reservoir for which the t = oo reduced density 
matrix of the model of Section 2 depends on p(0). 

7.2. The regime of an almost flat density of states of the reservoir. 

There is a case where the pa,a{E, oo) are equal to 1/2 and the off diagonal entries vanish: 
thus corresponding in some way to the Gibbs distribution ()3.21|) at infinite temperature 
(3 = 0. This occurs for an "almost fiat" density of states z/q of the reservoir. More precisely, 
we assume that has the form (cf ()H.19|l ) 

MX) = -if (-\ , (7.3) 



where is a non-negative function, having unit integral, and a is the biggest parameter 
having the dimension of energy of the problem. We assume in addition that the parameter 



^-^ = .M^y:=A (7.4) 
a 

is fixed, i.e. that the coupling constant is of the order y/a . Hence, the regime corresponds 
formally to the following limiting procedure: 

a — > oo, — > oo, A = 7ry9(0)f^/a is fixed. (7.5) 

By using the relation 

dx 



7ria{z), a{z) = siga{^z), (7.6) 

X — z 

we obtain that in this regime 

v\^{z) = Ma{z), r^iE, z) = ^ ^ _\ _ j^^^^^y (7-7) 
Thus, according to (jSiHl) and (Hl!^ . we have (cf (gUHl)) 

^^^^^^) = ^TV—-^^^vr—^V t;V,(A) = A/vr. (7.8) 
Ti\\E + sa — A)"' + 

Now, computing respective integrals in ()4.9p (see Appendix 2), we find that 

p„,.(£.i) = ^ + a ""^'°';''-'"' e-". (7.9) 

where 

r = 4A = 47rt;Vo(0). (7.10) 
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Formula (|7.1(J|) for the rate of the exponentially fast convergence to the infinite time 
(equilibrium) limit (jT.lip is an analog of the Fermi Golden Rule in our context. 

These formulas can be obtained as a particular case of ()5.7p . and ()5.8|) . corresponding 
to the replacement of all arguments in i^o by zeros, which is natural in view of ()7.3|) and 
the spirit of the approximation, according to which a is the biggest parameters with the 
dimension of energy, i.e., a » E,sa, hence I'oiE ± 2sa) ~ i^o{0). Note, however, that 
(|5.7|) . and (j5.8p correspond to the small coupling case (jS.lj) . while (j7.9p . and (j7.1Up are 
obtained assuming, at least formally, that v ^ oo (see fl7.5|) ) It can also be shown that 
in this asymptotic regime the off-diagonal entry p^^^{E,t) of the reduced density matrix 
decay exponentially fast as t ^ oo with the rate 2{A — \J — s^), if < s < A, and with 
the rate 2A, \i s> A. 

It follows then from (j7.9j) that in this approximation the t ^ oo limit of the reduced 
density matrix is 

p(E, oo) = f °, V (7.11) 



1/2 

i.e. it has the "infinite temperature" (p(O)-independent) Gibbs form of ()3.21|1 . 
Likewise, we have for the transition probabilities from ()4.8|1 . and ()7.9|) : 

p„,,(E,t) = i + ^e-^*, (7.12) 



p(i5,oo)= ( \'l \'l ). (7.13) 



in particular, 

1/2 1/2 

Formula ()7.9|) is derived in Appendix 2. The derivation makes explicit a mathematical 
mechanism of the exponential decay of the second term (exponentially fast convergence 
to equilibrium). Namely, to compute the double integral in ()4.1H) with the functions 
ra{z) and ra{E,z), given by ()7.7|1 . we make analytic continuation of ra{Xi — iO) to the 
upper half-plane and of rQ,( A2 + ^0) to the lower half-plane, which are the "non-physical" 
sheets for these functions. This procedure can be viewed as a simple version of analytic 
continuation of the resolvent of a self-adjoint operator via the cut, determined by its 
spectrum, i.e. as a "toy" case of analysis of resonances. In the recent rigorous studies 
PI in] of the return to equilibrium this analysis was carried out by using the techniques of 
the complex dilatation, applied to the Liouville operator of the composite system. 

Formula (j7.9|) implies that the pa,a{E, t) satisfy the system of two differential equations 

^Pa,a = ^(Pa,a " P-a-a), a = ±. (7.14) 

Eq. ()7.14j) describes a two-state Markov process, known as the random telegraph signal 
or the dichotomy process (see e.g. [ID])- The process assumes alternatively two values 
on intervals, whose lengths are independent identically distributed random variables with 
the probability density Fe"'"^. The equation plays the role of a master equation for our 
model in the asymptotic regime ()7.3p - ()7.4|) of the almost flat density of states of the 
reservoir. It is the (3 = version of the equation 1)6.111) . 

7.3. Special cases of the generalized model. 
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We consider here certain asymptotic cases of the model of Section 6, defined by 
Namely, we assume that the density of states uq of the reservoir has the "almost flat" 
form (Q - (d), that the form-factor f{E, E') of is 

f(E,E')=w(^^], (7.15) 



and that, in addition to ()7.4|) . the inequality 

a > 6 (7.16) 
is valid. This leads to the following system of equations for the self-energies Aq, of ()6.4j) : 

^"(^-) = -4 r /"(-z")"; (7.17) 

J-oo x + E + sa-z - A_a{E, z) 

where A is defined in (j7.4j) . To compute the integral we take the convenient "Lorenzian" 
form of w 

w{x) = (7.18) 

and we obtain 

TV Ah 

^c.{E,z) = -- (7.19) 

E + sa — z — A_a[E, z) — iba[z) 

where a{z) is defined in ()7.6|) . Similar argument, applied to the operator Sa,s of ()6.5|] . 
shows that ()6.2|) has the form ()4.3|) in which v^Ta is replaced by from ()7.19|) . 

Before analyzing the reduced density matrix in this asymptotic regime, consider the 
case, where b ^ s, i.e. where the system dynamics is negligible. In this case the r.h.s of 
(17. 19|1 coincides with the r.h.s. of ()7.7|1 . thus the diagonal entries of the reduced density 
matrix are given by ()7.9|) . This result for the transition probability was obtained in ^3] 
by using the Dyson series in powers of the system-reservoir interaction for the evolution 
operator exp{—it -ffnL=o} keeping only those terms of the respective series that are 
relevant in the n ^ oo limit. They considered the regime defined by ()7.3|) . ()7.4|) . and 
used the inequalities a ^ b ^ s. 

On the other hand, in the regime (|7.16p but with the energy s of the small system 
being of the same order as the range of energies b, participating in the system-reservoir 
interaction, we obtain from ()7.19p 

By using (j7.2(jp and the above remark on the form of (|6.2|) . we obtain again formula (|7.9|) 
in which T of ()7.10|) is replaced by 

ri = rj^. (7.21) 

This result was obtained in |16j, by using a more sophisticated version of the perturbation 
theory, proposed in 
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Appendices 



1. Large energy tails of //^(A) 

We will prove here the asymptotic formula ()A.9|) . Recall that the density Ua is defined 
via its Stieltjes transform 

r.{z)=l'-^,Qz^O, (A.l) 

by formula ()3.8|1 . and by the system ()3.9|1 . Because of formula ()3.8|) it suffices to consider 
the case, where '^z > 0. By using the identity 



1 



A-C 



e-'^^-'^^'dt "^z > 0, 



and taking into account that, according to ()A.1|) . '^ra{z) ■ '^z > 0, '^z 0, we can write 
dnSl) for Qz>0 as 

POO 

r^{z) = % / e-'*(""-"-"'''-"("))/jWt^^, (A.2) 

where 

fj{z) = J e-''^uoiE)dE = exp l-ttcoJ - ^ j . (A.3) 



In writing the last equality we used the form ()3.19|) of the density of states of the reser- 
voir. Since Qra{E + iO) is non-negative for '^z > 0, formulas ()A.2|1 and ()A.3jl imply the 
inequality 

K{E + tO)\ < {n /2Ja^y^^ (A.4) 
The same formulas lead to the relation 

2 2 2 P 2 

Jo 

where 

^ = x/a^, X = e — Cq. (A. 6) 



Replace the integral in the r.h.s. of ()A.5|) by the sum of the two integrals Ji and I2 over 
the intervals [0, i^] and [i^, 00 + ^y+) respectively. The first integral can be written as 

Jo 

where ^ 

^i(^) = ^ — 7^=' 92{u) 



1 + 2vr^ 

In view of flA.4|) the leading contribution to Ji as J — > cxd is due to the lower integration 
limit u = 0. This yields the asymptotic formula 



where a; and ^ are defined in ()A.6|) . 
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The integral I2 can be transformed to the form 

2 2 /"OO 2 2 



The leading contribution to I2 is also due to the lower integration limit u = 0, hence we 
have 

2 2 /* 2 2 

/2(e) = .e-^e-^^^'^-'^^-^^+^o" / e-^dM(l + 0(J-i)) 



^2 2 



2Ja2' 



These formulas imply (cf ()A.4|) ) that ra{Je + iO) has the asymptotic form, coinciding with 
the r.h.s. of ()A.7|1 . and, as a result, we obtain that 

z/„( Je) = 7r-i53r«( Je + ^0) = ""'^^ ^j^^ ^~.{Je) + z/o( Je)e«^"(l + 0( J-^)), (A.8) 

where the first and the second terms in the r.h.s. are the contributions of 7r^^$5/i(Je + zO) 
and 7r~^53/2( Je + iO) respectively and z/q is given by (j3.19j) . By combining ()A.8|) for a = ± 
and ()A.6|) . we obtain 

z/q,( Je)/z/o(-'^e) ^ exp |(e — eo)sQ:/a^} , J ^ 00. (A. 9) 

This formula and (HnUll yield (jSTTj). 

Furthermore, formulas ()4.13|) . ()4.4|) . and ()A.8|) imply that for any continuous function 
(f) we have 

/oo 
0(e)i/„(Je, Je') Jrfe' = 0(e), (A.IO) 
-00 

i.e. that Vai.J^i Je')J tends to 6{e — e') as a generalized function. 
^. Derivation of formula (|7.^| ) 

We outline here a derivation of formula (|7.9|) for the mean reduced density matrix of the 
small system in the regime (j7.3p - (j7.4|) of almost flat density of states. By using formulas 
fj4.10|) and ()7.7p for ra{z) and ra{E,z), we can write ()7.9|) as 

PaAE,t) = -J^J Ai - Ia "^^' J e"'*^'[(^i-^2-2iv4)p,,jO)-2iAp_,,_,(0)] 

(Ai - A2 - 2tA) 

(Ai-A2)(A2 + M)(Ai-A2-42A) 

According to ()4.10|) . the integrand of the internal integral in A2 is the limiting value of a 
function that is analytic in the upper half-plane and has a cut on the real axis. We see, 
however, that the function can be continued in the lower half-plane and has there simple 
poles at A2 = Ai — iO, —iA, Ai — AiA. Computing the integral by residues, we obtain 
after certain transformations: 

1 r Jt^i^ - 2^^)Pa,a(0) - 2iAp_«,_«(0) 



2m - AiA) 
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where the integration path encircles ^ = from below. The integrand can be continued 
in the upper half-plane, and has there the simple poles at ^ = + iO, HA. Computing 
the integrals by residues, we obtain ()7.9|) . 

Our argument above is a toy version of a techniques of complex dilatation, used in jHlC] 
to find the resonances (poles on the second sheet) of the Liouville operator, responsible for 
the exponential convergence of the density matrix of the composite system to the Gibbs 
distribution. Hence, the mathematical mechanism of the exponential convergence in our 
model is similar to that of the spin-boson and related models, studied in [HlC]. 
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